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Abstract 

We present all-order results for the highest three large-jc logarithms of the splitting functions P qg 
and P gq and of the coefficient functions C^ q , C2, g and C^ g for structure functions in Higgs- and 
gauge-boson exchange DIS in massless perturbative QCD. The corresponding coefficients have 
been derived by studying the unfactorized partonic structure functions in dimensional regulariza- 
tion independently in terms of their iterative structure and in terms of the constraints imposed by the 
functional forms of the real- and virtual-emission contributions together with their Kinoshita-Lee- 
Nauenberg cancellations required by the mass-factorization theorem. The numerical resummation 
corrections are small for the splitting functions, but partly very large for the coefficient functions. 
The highest two (three for C^g) logarithms can be resummed in a closed form in terms of new 
special functions recently introduced in the context of the resummation of the leading logarithms. 



1 Introduction 



The splitting functions governing the scale dependence of the parton densities of hadrons and the 
coefficient functions for inclusive deep-inelastic scattering (DIS) are benchmark quantities of per- 
turbative QCD H). At this point these are the only quantities depending on a dimensionless scaling 
variable, the parton momentum fraction and the Bjorken variable (both usually denoted by x), for 
which third-order corrections in the strong coupling constant a s are fully known. The correspond- 
ing three-loop calculations started with DIS sum rules [J2J-3J and proceeded via a series of low 
integer moments of the splitting functions and the coefficient functions for the most important 
structure functions [EH3 to the corresponding complete calculations of Refs. II8HT21I. 

Such higher-order calculations do not only improve the numerical accuracy of the predictions 
of perturbative QCD but also help to uncover general structures, for example in the soft-gluon limit 
x — > 1. Writing the expansion of the splitting functions in the MS scheme as 

00 

P ik (x,a s ) = £a s M+1 4 n) (x) with a 8 = -i, (1.1) 

n=0 471 

the diagonal (quark-quark and gluon-gluon) splitting functions take the form 

Pt%) = 4° (I-*); 1 +*£°8(1-*) +C«ln(l-x) + 0(1) (1.2) 

(n) (n) 

with the n-loop quark and gluon cusp anomalous dimensions related by A g /A q = Ca/Cf at 
n < 3, where Ca and Cf are the usual SU(N) colour factors with Ca = 3 and Cf = 4/3 in QCD lfT3l . 
It was not known before Ref. Ifl4l . inspired partly by observations made for the three-loop results 
in Refs. (HUH, that the third term in Eq. (11.21 ) is linear in ln(l— x) at all orders n, and that its 
coefficients cj~ (with C% = 0) are simple functions of lower-order cusp anomalous dimensions. 

The form of the off-diagonal (quark-gluon and gluon-quark) splitting functions, on the other 
hand, is not stable under higher-order corrections but shows a double-logarithmic enhancement, 

2n-l 

41« = I 4*' } ln 2 - f (l-x) + 0(1) . (1.3) 

£=0 

The terms with £ = form the leading-logarithmic (LL) large-* approximation, those with £ = 1 
the next-to-leading-logarithmic (NLL) approximation etc. Recently an all-order resummation of 
the former contributions to Eq. (11.3b has been presented [ 15]. A main purpose of this article is to 
extend those results to the £ = 2 next-to-next-to-leading logarithmic (NNLL) terms. 

The dominant ( 1 — x)^ 1 large-* contributions to the quark coefficient functions for gauge-boson 
exchange structure functions in DIS such as F2 and F3 (and the gluon coefficient function for 
Higgs-exchange structure function F* in the heavy-top limit lfT6l[T7l0 also show a double log- 
arithmic enhancement. These terms are resummed to all orders by the soft-gluon exponentia- 
tion [fT8] - [2"5]| which presently fixes the coefficients of the highest six logarithms analytically, and 
the seventh term for all numerical purposes since the effect to the presently unknown four-loop 
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cusp anomalous dimension can be neglected in this context I125W271 . A resummation of the highest 
three (l—x) logarithms has been inferred in Ref. [|28l from the properties of the corresponding 
flavour non-singlet physical evolution kernels. While those subleading contributions to the 'diago- 
nal' 0(a® ) coefficient functions are not the main topic of this article, we will be able to verify the 
DIS part of those results and fix the only missing coefficient for the fourth (N 3 L) logarithms. 

The 'off-diagonal' 0(a\) coefficient functions, such as 

oo 

C at k(x,a s ) = £a" c i"*M for a, k = 2, g or 4>, q , (1.4) 

n=\ 

receive a double-logarithmic higher-order enhancement as x — > 1 as well, 

2n—2 

£<f In^-'Cl-oO + oU). (1.5) 

£=0 

Also here the I = LL coefficients have been determined at all orders n in Ref. Ifl5ll . and also here 
we will extend those results by deriving the corresponding NLL and NNLL results. 

Finally we will also address the coefficient functions C^.q and Q, )g for the longitudinal structure 
function Ft- These quantities have a perturbative expansion of the form (|1.4I) . and are given by 

= D ( l,? ln 2n - 2 ~ e (l-x) + 0(1)) , * = q,g, (1.6) 

V £=0 J 

at large x, i.e., they are suppressed by one power of (l—x) and In (l—x) with respect to their 

counterparts for the structure function F2. The coefficients D^q for £ = 0, 1,2 have been obtained 
already in Ref. 11291 , again from physical-kernel considerations. For the gluon coefficient function, 
however, only the LL coefficients ^ have been determined completely in that article. Below 
we will verify those results and extend also the resummation of Ct.g to the NNLL terms. 

The remainder of this article is organized as follows: In Section 2 we derive all-order expres- 
sions for the (dimensionally regulated) Mellin-/V space transition functions as far as required for 
the mass-factorization of the structure functions at the level of the dominant and sub-dominant 
contributions discussed above. The NNLL resummations of the unfactorized partonic stucture 
functions T a ^ are then derived in Section 3 for the cases (|1.4I) . where to NLL accuract we employ 
two different methods to drive the same results, and in Section 4 for Ft- The partly rather lengthy 
results of these three sections are then combined, and in Sections 5 and 6 we present and dis- 
cuss the respective resummed expressions for the moments of the off-diagonal splitting functions 
- which receive rather small resummation corrections at relevant values of iV - and the coefficient 
functions (11.41) and Ct, g for which these corrections are (very) large. We summarize our findings 
in Section 7 where we also present a brief outlook to future extensions and applications of some 
of the results. Closed expressions have not been found so far for the third logarithms in Eqs. (11.31) 
and (11.51) . Numerical and symbolic tables of NNLL coefficients to high orders are therefore finally 
presented in Appendix A for the splitting functions and Appendix B for the coefficient functions. 
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2 Large-* / large-N mass factorization to all orders 



The main part of our calculations is performed after transformation to Mellin-Af space, 



f(N) 



f l dxx N - l f{x) or /(JV) = [ l dx(x N - l -l)f(x)+, (2.1) 
Jo Jo 

where the ubiquitous x-space Mellin convolutions are reduced to simple products. To the accuracy 
required below, the relations between the large-* logarithms and their moment-space counterparts 
are given by 

{ - 1)k ( ^ k ~i-x X) ) + = l(^-m k + \k(k-i% 2 [suN)] k - 2 

+ i^-l)(fc-2)C 3 [Si-(^V)]^ 3 + 0(V(AO]^ 4 )) , 

(-l) k \n k (l-x) = ^(\n k N + h(k-l)^ 2 ln k - 2 N + lk(k-l)(k-2)^ 3 \n k - 3 N 

N \ 2 6 



+ 0\hi k - 4 Nj) , (2.2) 
(-l)*(l-x)ln*(l-x) U ^ ( ln ^ " k\n k - l N + ^k{k-\% 2 \n k - 2 N + o(ln fe - 3 /v)) 

withSi_(A0 = InA 7 - 1/(2N) + 0(l/N 2 ) and N = Ne^, i.e., ln/V = lnAf + y e with y e ~ 0.577216. 
Here = indicates that the right-hand- side is the Mellin transform (12.11) of the previous expression. 

The primary objects of our resummations are the dimensionally regulated unfactorized par- 
tonic structure functions or forward Compton amplitudes T a ^ for the combinations of a and k of 
Eqs. (11.41) and (|1.6I) . For brevity suppressing all functional dependences on N, a s and the dimen- 
sional offset £ with D = 4 — 2e, these quantities can be factorized as 



Ca,iZik ■ (2.3) 



Here the process-dependent D-dimensional coefficient functions (Wilson coefficients) C a i include 
contributions with all non-negative powers of £. The universal transition functions (renormalization 
constants) collect all negative powers of £ and are related to the splitting functions in Eq. (11.11) 
(or the anomalous dimensions y) by 



dZ 
d~\nQ : 



-y = P = -hzt^z- 1 . (2.4) 



Here and below we identify, as already in the introduction, the renormalization and factorization 
scale with the physical hard scale Q 2 without loss of information. Using the D-dimensional evo- 
lution of the coupling, 

d % = -ea s + p> s ) (2.5) 



din g z 

where P(a s ) denotes the usual four-dimensional beta function of QCD, P(fl s ) = — Po^s + ■ • • w i m 
Po = 11/3 Ca~ 2/3 rip Eq. (12.41 ) can be solved for Z order by order in a s . 
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The general elements of Z'"J become extremely lengthy at very high powers n of a s . Here, 
however, we are interested only in the LL, NLL and NNLL contributions at order TV -1 for Z qg and 
Z gq and N° for Z qq and Z gg (required for Eq. (12.31) also in the off-diagonal cases). Consequently 
there can be at most one off-diagonal N~ l factor per term. Moreover pj^ in Eq. (|l.ll) can enter 

only (once) at NNLL accuracy, and higher-order coefficients P^ in Eq. (|l.ll) are not at all 
relevant at this level. Finally Po in Eq. (|2.5I) and P^ only contribute from the NLL and NNLL 
terms, respectively, and P„>i would enter only at the next logarithmic accuracy. All this can be 
easily read off already from the well-known third-order expression for Z, 



Z - '+^ l0,+ <2^ ,0, - P °K , + ^ <1 
+ 4(^(Y ,0) -Po)(Y ,0) -2P„)y<°> 

>-2Po)t <1) + (y (1) -P.)2 Y W 



+ 



(2.6) 



together with Eqs. (fL2l) . (TOT) and (TO) above. 



The terms that do contribute to the a ^ b off-diagonal entries of Z at the present accuracy can 
be grouped as follows: 



7 (k) _ (k) 

^ab — ^ab 



+ ^ab 



P< 



+ z (k) 
+ ^ab 



+ Z (k) 



7(1)7(1; 



, z (k) 



7(1) Y m 



(2.7) 



The first term on the right-hand- side collects all contributions with at most one higher-order anoma- 
lous dimension y^- 1 ) but no contribution from the beta function. It starts at k = 1 and reads 



r(k) 
^ab 



1 



kl£ k 




y^a) yil (lib 



+ e£i( t -;-2)( t -;- 1 )( Y «?)^r^ , (YS 



k-j-3 



+£ 'fi (t -i-2)(Hi + i)W') H " 7 i;»,« (~<°»v 



(2.8) 



The first line includes, for i = 0, the LL expression used in Ref. |fT5l . The contributions linear in 
Po contribute from k = 2 and NLL accuracy and are given by 



7 {k) 
^ab 



Po 



2 kle k 



i=0 



7=0 



(2.9) 



while the corresponding NNLL Pq term in Eq. (12.71) for k > 3 is 



^ab 



Po 2 



i^ll e ' *X; , ^^[A-U-l)(A-2)(3A-I)-6/(/ + ./+I)A-(A-l; 



24 tie* ^ y f j! 



+ i(3/+l)(i + j + l)(/ + 7 + 2)](y < 



(o) 

Cttl 



(2.10) 



Finally we distinguish NNLL contributions with y( 1 )y( 1 ), where one of the factors is the off- 
diagonal 0(N~ l ) entry, and contributions with yWy^ 1 ), where the latter has to be off-diagonal 
at the present level of accuracy (recall that every term includes one off-diagonal anomalous dimen- 
sion J a ^b). The former terms contribute from order and are given by 



^ab 



' ' *" 7=0 j=0 



(y^'yS (i£»)'iS> (r£ , ) i " i " y " 4 +(^ ) )'^ ) W^' 

and the corresponding final contribution to Eq. (12.71 ) at k > 5 reads 



1 *-3 i+l fk-^k-i^t-l^ _ ^ ( k _i_j_fy 



Y 0) Y w We* ^ 2 | f t- ; 0) ; (fe-i-7-£-3)! 

k-i-j-£-3 



Jaa^j Jaa (jaa^j jjfi (ibb 



7 {k) 



k-£-3k-3-i-e 

The coefficients Eqs. (|2.8I) - (12.121) have been inferred by analyzing the respective first five to 
seven non-trivial orders k and then verified to 'all' orders using, as for a large part of our symbolic 
manipulations, the programs Form and TForm H30U3U . 

The corresponding result for the diagonal entries of the Z-matrices are much simpler due to 
Eq. (I1.2I ). Including also terms which contribute to the next-to-next-to-next-to-leading logarithmic 
(N 3 LL) terms suppressed by one power of 1 /N, the coefficients at order ocf are given by 

Laa ~ k\ v aa J 2{k-i)\ v aa ) Jaa 2 {k-2)\ v aa ) 

Po ^ ^ ( (0)\* ^ & ( (o)\*~ 3 (2) 

+ 24 Jk^yy. {3k ~ 1} v aa ) + 3(k-3)\ v aa ) Jaa 

12 (k-3)\ {3 b) V aa ) Jaa 3 (Jfc-3)! V m ) 

+ £i (^ , )'" , (^) J+ !S lt - i),3t - !) W 4 ^ 

Only the first four terms contribute to the iV NNLL expression entering the off-diagonal N~ l 
mass factorization (I2.3I ). Note that, unlike at order A^°, pj^ and are NLL and N 3 LL quantities 
at order A^ 1 due to cf } = and d"- 2) ^ in Eq. Q). 
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The unfactorized structure functions (12.31 ) are given by these results multiplied by 

oo oo 

C a ,i = 5 fly 8 iq + 8^8 /g + £ < £ z k Caf ] ■ (2-14) 

n=\ k=0 

The index y generically represents the gauge-boson exchange structure functions (except for F£). 
c i" ^ = c i"^ are me order coefficient functions in Eq. (|1.4I) for any combination of a and i. 
The quantities c^" - usually denoted by af], etc in fixed-order calculations - are enhanced 
by factors \n k N with respect to those four-dimensional coefficient functions. 

The calculation of T a to order af- n+1 and Z n ~ l (for F L : e"^ +1 ) provides the N"LO (leading- 
order, next-to-leading-order etc) renormalization-group improved fixed-order approximation to the 
structure functions F a . It is obvious from Eqs. (12.61) - (12.141) that a full N"LO result completely 
fixes the highest n+1 powers of 1/e to all orders in a s . An all-order resummation of the splitting 
functions and coefficient functions requires, at the logarithmic accuracy under consideration, an 
extension of these results to all powers of £. The flavour- singlet structure functions considered 
here are fully known at N 2 LO from Refs. ll8l- [T0l[T6l[T7l and the earlier coefficient-function cal- 
culations of Refs. [|32] - l35l . Hence a double-logarithmic resummation based on these results can 
be expected to predict up to the highest three logarithms at all higher orders, including the corre- 
sponding contributions to the three-loop coefficient functions for F2 and F§ exactly computed in 
Refs. IfTTlfm and the large-x predictions of the four-loop splitting functions in the latter article. 



3 Resummation of the unfactorized expressions for F2 and F§ 

In this section we derive all-order expressions for the leading ./V -1 contributions to the off-diagonal 
amplitudes or unfactorized structure functions Tz, g and T^ q at NNLL accuracy. We will apply 
two approaches: first an iteration of amplitudes generalizing the leading logarithmic results of 
Ref. [15], and then an apparently new and more rigorous treatment which makes use of only the 
D-dimensional structure of the unfactorized structure functions in the large-x limit and the KLN 
cancellations Il36ll37l between its real- and virtual-emission contributions. 

Both calculations require the corresponding expressions for the N° parts of Z qq , Z gg , C2, q and 
C^g which can be determined from the diagonal amplitudes T^a and 7^ g in the limit governed by 
the soft-gluon exponentiation. These quantities are given by 

T a , k = exp (a s I $ + al T® + a s 3 T® + . . .) (3.1) 

with 

oo 

T$ = I «f exp(nelnAO - vjf) . (3.2) 

C=-n-l 

For the quark case the coefficients entering the highest four logarithms at all orders in a s and £ read 

4V 2)=4 ' 4q _1)=3 ' 4^ = 7 ~ 4 ^ 4^ = 14-3^ 2 -8^ 3 , 

^=4, ^=6, V™ = 16 + 2&, ^ = 32-3C 2 -f C3, (3.3) 
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-(2,-3) _ T/ (2-3) _ o 
K 2,q ~ V 2,q ~ Po > 

p (2-2) _ /4 ?r \ r ^ 19 R u (2-2) _ (2,-2) 3 

< q _1) = (^-6C2 + 12C3)c F +(^-20C 3 )c A +(^-3C2)Po, 

= (^-12^ + 24^)c f -(^ + 26^)Q-(^ + ^ 2 )Po, (3.4) 

n(3-4) _ „(3,-4) _ 4 2 v (3-3) _ -(3,-3) R 2 

R 2,q ~ V 2,q ~ 9 Po > V 2,q ~ ^2,q + Po > 

<" 3) = ~Ci-fc F C A + ?C f p„ + (|-^fc) CA p„ + |pJ (3.5) 

fi <4,-5) = ^(4,-5, = lp3 > (36) 

where we have suppressed an obvious overall factor of Cf and expressed the dependence of the 
number of effectively massless flavours in terms of Po- To N 3 LL accuracy these results are 
converted to the renormalized coupling used elsewhere in this article via 

d 8 = (g_!) a 3 + | a 4 + .... (3.7) 
The corresponding gluonic coefficients are required only to NNLL accuracy here and read 

(|-4C2)c A + |po, 

2C 2 Ca, (3.8) 



and 





-2) 


= 4C A , 




1,-1) -(1,0) 

•,g -P°' \e 




-2) 


= 4C A , 




;i,-i) n (1,0) 
,g _ u ' \g 


d(2,- 
\g 


-3) 


<t>,g 


3 ) = 




r,(V 
\i 


-2) 


-(!- 


■2C2 




^,g 


-4) 


= Vif 3 '" 

<t>,g 


4 ) = 





C -C-jW' (3 - 9) 

(3.10) 

After the transformation to the renormalized coupling, Eqs. (13.11) and (13.21) for T^g need to be 
multiplied by the renormalization constant of G MV G^ V D81I39L 

1 - 2p\ ) e~ 1 tf 2 + 3p 2 e _2 a;? + ... . (3.11) 

Note that there is no new physics content in Eqs. (13.11 ) - (13.101 ) which represent the Mellin 
transforms of the decomposition of T2, q and T§ :g used in Refs. [|40ll4T1l and its obvious higher- 
order generalizations, cf. Ref. (421, recast in an exponential all-order form. Recall also, from the 
same references, that of the £ coefficients /jw -3 ), ^( 2 ^-5) ^(n,£-2n-l) re i evan t f or the £-th 
logarithm at order n only one combination is not fixed by lower-order information, and that the 
same holds for their virtual-correction counterparts V (M -3 ), V ( 2 ^~ 5 ) . .. V ( n ^ _2 »-l) 
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We are now ready to address the resummation of the off-diagonal amplitudes. Using the par- 
ticularly simply colour structure of the identical leading-logarithmic contributions to 7^ q/Cp and 
T2,g/rip their LL resummation has been inferred in Ref. lfl~5il to which the reader is referred for a 
detailed discussion. The result can be written as 

r (») ll 1 T (i) 'v ( \ ~ X r^ r( n_i_1 ) ex \i\ 

L a,k = n 1 a,k L\ i ) >,g 7 2,q ^ A£ > 

where, of course, only the respective first terms of 

T$ = -^exp( e in^V)(2-£-(3-2C2)e 2 + ...), 

r zg = -^exp(£ln^V)(2-2£-(6-2^ 2 )£ 2 + ...) (3.13) 

are required. The corresponding expressions for T^q and 7^ g can be read off from Eqs. (13.11) . (13.21) . 
(1331) and (gSJ) above. 

Eq. (13.121) can be generalized to next-to-leading logarithmic accuracy, where running-coupling 
effects enter for the first time, using the natural ansatz 

7$ N = L k<]>( E/^i.BjrWrf- 1 ) - £ E%^(»,0r«^- 2 )} (3.14) 

where f a ,k{ n , i, £) are linear functions of £. It is not a priori clear that such a simple ansatz is com- 
patible with the infinite number of constraints provided by the NLL coefficients of the four highest 
powers of l/£ at all orders n in a s which are provided by the results of Refs. JHHMIIMITl mcm ding 
the large-x results for the four-loop splitting functions Pq g and Pg q in Eq. (|l.ll) . However, all these 
constraints can indeed be fulfilled, and the resulting coefficients are given by 

f^ q (n,i,e) = ( n ; ) 

fi, g (n,i,£) = ("y 1 ) 

and 

S^q(M = g 2 , g M = [ill ) ( 3 - 16 ) 

with Bkj = 1 for k > j and 0^- = else. Eqs. (|3.14l) - (|3.16l) together with their diagonal counter- 
parts (13.11) . keeping the £ 2 and £ 1 contributions to T^) and the £ 3 terms of in Eq. (|3.2I) . 
facilitate the extension of the all-order mass factorization to the next-to-leading logarithms. 

An extension of Eq. (13.141) to the the third logarithms can be expected to become much more 
cumbersome, requiring at least £ 2 corrections to Eqs. (13.151) . £ corrections to Eqs. (13.161) and a new 

contribution, but presumably also terms respectively involving tS 1 ] or TPlj . Instead of pursuing 
this approach, we now switch to our second method for the resummation of T^g and T^q. 



1 +£ 



A(;+l)(„-0e zl -^(l-n8 /o ) 



1 + £ 



8C A 



i(n-i-3) +-(3j + l-n9„28/„-i) 



(3.15) 
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For this purpose we consider the calculation of T a £ via suitably projected gauge-boson parton 
cross sections as performed at two loops in Refs. Il32] - [34l . The maximal (2 — >■ n + 1 particles) 
phase space for these processes at order a" is H43[I44| 

(\-xy l ~ n£ f d(3n-2 other variables) f(x, ...) , (3.17) 
Jo 

where one trivial azimuthal integration has not been counted. The integrals for the n-th order 
purely real (tree graph) contributions T^j do not lead to any further factors (1— x)~ £ , hence their 
expansion around x = 1 can be written as 

T ( " )R = (l-x)- 1 -"' V (l- X f^—{ R ^ + £i?WNLL E 2 R (n)NML \ 
a,j v i V > p2n-\ 1 a,],% a,],% a,],% I v ' 

The mixed contributions (2 — >• r + 1 particles with n — r > 1 loops) include up to n — r additional 
factors of (1 — x)~ e from the loop integrals on top the phase-space factor, leading to 



e=r ^=0 



• iu WLL + £M (n)NLL + £ 2 M W ™ L + 1 (3 19) 

Finally the diagonal cases, where terms with \ = are present in Eqs. (13.181) and (13.191 ), also 
receive purely virtual contributions T a y V given 
to an amply sufficient accuracy 1140114 l[|45l - |4~8l 



receive purely virtual contributions T a y V given by the y*gg and form factors which are known 



r fl ( f = 5(i-x)^{yjf L + £yjf NLL + £ 2 v ( 'f NL + ...). (3.20) 

Note that for ^ = the (1— x) factors in Eqs. (|3.18l) and (13.191) are D-dimensional +-distributions 
which include a factor £ _1 8(1 — x) after expansion in £. 

(n) — 

The partonic cross sections T^ k J at order a" do not include any poles higher than £ , i.e., 
T$ = + r fl f + (8 flY 6 ( - q + 8fl+ 8, g ) rjf 



L{r {n 

£ n \ a,j 



Y> + E T<$ 1 +BT<3 2 + ...\ (3.21) 



(recall Eq. (12.141 ) concerning the index y). Hence there are n — 1 — m KLN relations between 
the n N m LL coefficients in Eqs. (|3.18l) and (13.191) . Since a full N fc LO calculation provides, as 
discussed above, the first k+l non-trivial powers of£~ ! ,it leads to a total of n — m + k relations. 
Consequently the coefficients up to the N^LL terms are fixed (and those for m < k over-constrained) 
in terms of the N^LO results by nothing but the above D-dimensional structure and the mass- 
factorization formula (|2.3I) which is guaranteed for structure functions in DIS by the operator- 
product expansion ||49l , see also, e.g., Refs. [15011511 . 
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Note that this resummation is far less predictive than the soft-gluon exponentiation 11181 - 1221 of 
the (1— x)~ /N° terms of C2, q and C^g which involves an additional factorization in the threshold 
limit. Therefore, as mentioned below Eq. (13.1 II) . only one of the n N m LL coefficients for £, = is 
actually independent. Hence a N"LO calculation in this case implies a N"LL exponentiation which 
fixes the highest 2n+l logarithms at all orders in oc s . 

Now we switch back to Mellin space and apply the above results to and T^q. The n-th 
order contributions to these quantities can be written as 



T$W = (A^+BB^+^cSf + ...)exp( e („-/)lnAr). (3.22) 

(=0 



As discussed below Eq. (|3.21l) . only one of the LL coefficients A I' is independent for each value 
of n and a, k. Choosing i = for that, the KLN constraints on the other coefficients read 

<' t ' , = (-i) i ("7 1 k:f ) - 023) 



The i = coefficients are found to be 

1 , («,0) _ 1 A (n,0) _ n 2n-\ 1 \> r £ r n-£-l CX1A\ 

At NLL level two coefficients in the sum in Eq. (13.221) are not fixed by the KLN cancellations. 
In terms of our choice i = and i = 1 the remaining coefficients are given by 



B (n,i + D = { _ iy 



a ' k \ i + l a ' k 



(3.25) 



The all-order results for the respective two coefficients determined by the NLO results are 
fl(j°> = -^n/{ n fc^- £ - l (nn 2 + 39n + 22-Ue) 



+ 36C^ 1 (n + l)0„2 + C£ _1 (lln 2 + 15n + 22) 

n-2 



n-l \ 

-2n f Y,Cfc£- e - 2 (n 2 -3n + 2(£+l)) \, (3.26) 
e=o ' > 

B t% = ^«/{'E c / c A^ 1 ( lln3 + 50n2 -( 9 + 14 (^- 1 )) n - 66 + 6 ^) 

on. I f-i 

+ 36C^ l (n + 2)(n-l) + C n k ~ x (l In 3 + 26n z + 29n - 66) 
-2n f Y,Cfc£- e - 2 (n 3 -2n 2 + (7 + 2£)n-6(£+\)) \ (3.27) 



10 



and 



B 



(n,0) 



B 



(n,l) 



^— C F I £ C|C^ _1 (1 In 2 + 33n - 70 + 14£) 

+ 12C;- 1 (9n-7) - 2n f Y,C^- e - 2 (n 2 + 3n-2(£+l))\ , 

-2 (- n-2 

— C/7 <^ Y< C F C A^ X (H« 3 + 44n 2 + (14£-81)n + 26-22£) 



(3.28) 



+ 12C;- 1 (9n 2 -13n + 4) 



«-2 

-2n / £c/c;-^ 2 (n 3 +4n 2 -(3 + 2£)n-2(£+l)) \ . (3.29) 

Similar to Eqs. (12.81) - (12.121) . we have first inferred these results by analyzing a couple of or- 
ders (enough to over-constrain the numerator polynomials) and then verified them to 'all' orders. 
Eqs. (I3T221) - (13T291) lead to the same results as Eqs. (13T3T) - (13T61) above. 

At NNLL accuracy, finally, all but three coefficients (chosen as i = 0, 1 , 2) in the sum in 
Eq. (13.221) are specified by the vanishing of poles higher that £~" in T 2 g and 71^ • These co- 
efficients can be written as 



_,(n,/+2) 



n — 3 
i 



C (w , 2) + r 



n-2 
t + 1 



«-i 

i + 2 



C 



(»,o) 

a,k 



(3.30) 



The general expressions for the three independent coefficient at rather lengthy. However, especially 
since we will not be able to express the NNLL mass-factorized results in a closed form, they are 
presented here nevertheless in order to assist future research by others. The coefficients for T^g are 



C 



jc£ _1 (76n 2 + 60n + 8) - C^~% (30/7 2 + 30n - 12) 



2 -g " __ 6Z" 

+ 300C F 8„2 - 48C^- 1 C 2 ne„ 2 + PoC£~ 2 (3n 3 -4n 2 - 9n+ 10) 0„ 2 

+ C;" 1 (^n 2 + ^n+lll)e„3 + C F C^ 2 (28n 2 + 132/7 -68)0 



V 2 " ' 2 

£=1 

n-2 



[3 4 

[32" - 


29 3 
48" 


+ ( 




r 9 3 / 


9 

4 ? + 


17 

~2 





41 \ 2 



n3 

49 \ 1 



J n U £+ 48r + 8 (£ +3£+2) J 0n3 



6 n- 



31 



'+10 



e„ 3 



- £c;c;^^ 2 [6n 2 + 54n-12(£+l)]0 



n3 



£=l 
1=2 



58n 2 -(54£- 156)n + ^ 



179 



' + 8 



?n4 



(3.31) 



11 



C 2 ( " g 1} = ^L-^jc^- 1 (76n 3 - 24n 2 - 52n) - C^ 2 (30n 3 - 12n 2 - 6n - 12) 

+ 456C F 8„ 2 - C^^ 2 (48n 2 -72n-24)e„2 + (3 C^ 2 (3n 4 -n 3 - 21n 2 + 49n - 30) 0„ 2 

(27 1 89 \ 

—n 3 + 1 1 In 2 + —n - 225 j B n3 + C F C n A 2 (28n 3 + 96n 2 - 152n + 144) n3 

+ (^ 2+4<+ ^)"-f (<2+M+2) ; 



-A 



3 5 31^ 4 , f 3 13 \ 3 
8 + 96/ 



96 > 



n-2 

=i 



+ 1 foc/cr <~ 2 

21 - 33 



9«3 



» 4 -(?^)» 3 -(^-f)" 2 -(? 2 -^-«)» 

e„ 3 



58n 3 - (54£ - 90)n 2 + ( y £ 2 + ^£ - 154 )n - ^-£ 2 + 



37 

y 



n-2 

~ E C f C a" £ "^2 [6n 3 + 36n 2 -(12£ + 78)n + 12(£-l)]0„ 3 |> , 
e=i 1 



81 
2~ 



369 



9;j4 



(3.32) 



in— 2 



<^ C n A x (38n 4 - 92n 3 + 6n 2 + 56n - 8) 



>,2) _ 

- C2~ 1 ^2 ( 1 5n 4 - 42n 3 + 39n 2 - 48n + 36) - Cf 1 ^ 2 (24n 3 - 96n 2 + 72n + 48) 0„ 2 



n-2 



+ PoC 
147 



1 



47 



153 



-n 3 - -/T - — n J + — n z - 104/z + 50 



L2 



9«2 + Cp 1 



27 4 111 3 
^ n+ -2~ n 



2 771 

n 2 n + 366 

4 2 



0„ 3 + C F C^- 2 Un 4 + I6n 3 - UOn 2 + 2Un - 22o] n3 



B-3 

+ 2^ PO C F C A 

£=0 

1 



V 13 



5 / 3 n 31 \ 
« 5 + (^ — £ 



/ V8 192/ 



64 192 V 16 64/ V8 192 

+ EP0C/C 
9 



r£/~in— I— 2 
-A 



0«3 



(^ + 31 < -MZ)„^(lH < 2 + ^_ 101 )„_H£ < 2 + 65 <+50 
V 8 2 / V 8 8 / 4 4 



0«3 



n-2 



- E C^C2~ £ ~ 1 ^2 [3n 4 + 6n 3 - (6£ + 69)n 2 + ( 1 U + 12)n - 12(£ - 3)] n3 

n-2 



+ Ec/cr £ - 1 

t=2 



29n 4 -(27£+17)n 3 +(y£ 2 ' 361 • '^- 2 



— ^-162 )/r 
4 / 



-(^£ 2 -^£-164)n + 81£ 2 -293£-8] n4 | 



(3.33) 
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The corresponding results for T^ q read 



-(»,o) 

'4>,q 



in— 2 



6n\ Cf { C ^ 1 (^" 2 + ^" 21 ) + c r^2(48n 2 -144n + 48) 
- C n f 2 C A ( y« 2 - y « + 7) e„ 2 - C;- 2 C A ^2 (54« 2 + 6n - 60) 0„ 2 



4n 2 + 12n + --£ 2 



B-3 

0„ 3 - £ C/C A "-^^ 2 [6n 2 + 42n + \2t - 36] w3 



n-3 



+ £ PoC/C 

£=0 
n-3 



t n n-l-2 



£ + 2)n 2 + (^-£ + 6)n 



\ 4 
25 , i 



4 4 
1 



0„3 



(3.34) 



^(n,l) 

'0,q 



+ iPoC/C^- 3 [±n 4 + ^n 3 - ±(121 - 13)n 2 - -L(54£+ 121)n 
+ ^ 2 + 3£ + 2)]0 n3 } , 

Cf { C ^ 1 (^" 3-9 " 2 ~^ n + 3 ) + c r^2(48n 3 -192n 2 + 216n-72) 

- C^ 2 C A ( y n 3 - 34n 2 + - 65^ 0„ 2 - C^ 2 C A ^2 (54n 3 - 60n 2 - 78n + 84) 0„ 2 

2o /27 4 31 o 343 2 91 \ rt 
+ C;- 2 p ( y« 4 + y" 3 - — « 2 + y" - 32 j 0„2 

„ 43 i 

0n3 



in— 2 



6n! 



+ E Clef'- 1 [ 4 " 3 + ( I * + f ' - 28)* + - y € - 32 

n-3 

0n3 



£ d F C n A l ~ 1 C2 6n 3 + 24n 2 + ( 1 2t - 90)n - 1 2t + 60 
=0 L 

n-3 



+ £ PoC/c;-^ 2 [ ( \i + 2) n 3 + ( ^ - 1 2) n 2 - ( ^ 2 + It + 26) n 



/45 



■18£ 2 -22£-4 



0n3 



+ ^n 4 - t^(36£ - 37)n 3 - ^(216£ + 533)n 2 



1 

96' 



+ EPo 2 C/C A - / - 3 

£=0 

+ ^(^ 2 -6£ + ±)n + i(23£ 2 + 69£ + 46) 



1 

96 



1 



0n3 



(3.35) 



and 



C, 



(n,2) 

<t>,q 



m-2 



6n! 



. ,/ 183 4 111 3 531 2 201 

" w — n 4 n 5 n 2 A n + 42 

4 2 4 2 

13 



+ C n f l 1,2 (24n 4 - 144n 3 + 3 12n 2 - 288n + 96) - C n f 2 C A ( — n A - 33n 



13 



+—n 2 - —n + 42^ n2 - C^ 2 C A ^ 2 (21n 4 - 90n 3 + 21 n 2 + \AAn — 108) 0„ 2 

^n-2o ( 21 5 71 4 819 3 1413 2 M 93 \ 

+ C; 2 p ( yn 5 + — n 4 - —n 3 + —n 2 -92n + — j Q n2 



2n 4 -Sn 3 +( 2 ^£ 2 + ^-£-l0)n 2 



139 



n-21£ 2 + 3£ + 12 



Jn3 



M-3 

- £ Cfc%- £ -% [3n 4 + (6£ - 63)n 2 -(W- 144)n + 12£- 84] 0„ 3 

(|, +1 y + (^_ 16 y_(| / » + ^_ 7 y 



£=0 
f=0 



l r n-l-2 



135 „ 9 241 



8 



8 



36 )n + 



117 



165 



+ 12 



0«3 



n-3 



2 n t n n-l-?> 



+ EPo 2 c|c; 

1 -,2 



3 6 95 

— n 6 H . 

L 64 192 



n 5 - — ( \2l + 23)n 4 - — (2SS£ + 901)n 
64 192 



< 2 + ^(135£ 2 + 561 £ + 329)n - \{3U 2 + \05£ + 70)1 0„ 3 j ■ 
48 8 J J 



(3.36) 



4 Resummation of the unfactorized expressions for F L 



Up to one additional power of £ the moments of the unfactorized longitudinal structure functions 
Ti A and Ti g are built up from D-dimensional exponentials in the same way as T2. g and 7* „, 



1 



n-l 



N l+S kg£ 2n-2 ^ 



I M& + e<T + e 2 C^ !j + • • • ) exp(£ (n - i) InN) . (4.1) 



,{n,i) , c 2 r (n,i) 



-2n+l 



In these cases also the £T n poles vanish at order a" which compensates the absence of £ 
contributions. Consequently Eqs. (13.231 ), (13.251) and (13.301 ) are valid also for F^. As discussed in 
the introduction, our main objective for Fi is the resummation of C^g. As in the case of F2 the 
corresponding D-dimensional quark coefficient function, and thus Tl a , is also required for this. 

(n,0) 



Also here the LL coefficients A L ! are very simple and closely related, 



(n,0) 



)2/j 



CI 



(»,o) 



22n+l 



L >4 " [n-l)\~ r ' (n-l] 
The corresponding NLL contributions to Eq. (|4.1I) are given by 

4«-2 



(4.2) 



5 



;«,o) 

L,q 



Cr l {Mn-l)(n-2) + 16C A (n-l)(l-C 2 ) 
-4C F (9n-13) + 32C F ^ 2 («-1)} , 



(n-l)! 



(4.3) 
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= -^-^cr l {Mn 2 + 3n-6) + 16C A (»- 1)(1 -fe) 

-4C F (9n-16) + 32C F ^ 2 («-1)} (4.4) 



and 



5 



= ^riyT"/ C A _2 { 2 Po(«-l)(«-2)e„ 2 - %C F {2n-\)Q n2 + 6AC A n} 

' ! »/E C/CJ-* -1 8*3 , (4.5) 

Z=2 



(n-l)! 



B i M g 1} = -7^Tjy«/C A "- 2 { 2 Po(« 2 -« + 2)(n-l)0„ 2 - 8C F n(2n-3)0„ 2 

^n— 1 w— 1 

("-!)! ,2 



+ 64C A n(n-l)} + 2-^— - W/ £ C/C^- 1 (n-2)e B3 • (4.6) 



The NNLL terms for 7^ q read 



C ^ " 12(n-l)! C ^ l |jo U" + 12 n 8" 12 n + 7 J 

+ C A 2 (96n 2 - 288n + 192) (£ 2 - £ 3 ) + PoC A (l2n 3 - 84n + 72) ( 1 - £ 2 ) 

- PoCf [3(9 - 8^ 2 )« 3 - 32n 2 - (293 - 168£ 2 )n + 2(149 - 72£ 2 )] 

+ 6Cj [( 17 + 48^ 2 - 64^ 3 )n 2 - (463 + 1 12£ 2 - 544^ 3 )n + 510 + 32^ 2 - 480£ 3 ] 

- 8C F C A [(4 + 45^ 2 - 48^ 3 )« 2 - (134 + 1 1 1£ 2 - 276£ 3 )n + 130 + 66£ 2 - 228£ 3 ] } , (4.7) 

^ " 12(n-l)! CF l Po U + 24" 24" 24* + 12 * ™) 
+ Cl (96n 3 - 384n 2 + 480n - 192) (£ 2 - £ 3 ) 

+ p C A (I2n 4 + 12n 3 - 228n 2 + 420n — 216) (1 — 

- PoCf [3(9 - 8^ 2 )n 4 - 2(7 + 12^ 2 )n 3 - 19(29 - 24^ 2 )n 2 + 2(649 - 420^ 2 )n 

- 8(95 - 54£ 2 )] + 6Cj [(17 + 48£ 2 - 64^ 3 )n 3 - (534 + 1 12£ 2 - 608^ 3 )n 2 
+ ( 1227 + 80^ 2 - 1 120^ 3 )n - 710 - 16^ 2 + 576^ 3 ] 
-8C F C A [(4 + 45^ 2 -48^ 3 )n 3 -(152+144^ 2 -324^ 3 )n 2 

+(308 + 165^2 - 540^ 3 )n - 160 - 66^ 2 + 264^ 3 ] } (4.8) 

and 
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L ' q 12(n-l)! F l H0 V 16 48 16 48 4 12 / 

+ C; 2 (48n 4 - 288n 3 + 624n 2 - 576n + 192) (£ 2 - £ 3 ) 

+ PoCa (6n 5 + 6n 4 - 222n 3 + 7 14n 2 - 864n + 36o) ( 1 - £ 2 ) 

- PqCf ( | (9 - 8^ 2 )n 5 - 1(23 + 24^ 2 )n 4 -|(313- 296^ 2 )n 3 + |(1289 - 952^)n 2 

- 24(115 - 72^ 2 )n + 2(647 -360£ 2 )) 
+ 3Cj? [(17 + 48^ 2 - 64^ 3 )n 4 -2(311+ 80^ 2 - 368^ 3 )« 3 + (2675 + 144^ 2 - 2432^ 3 )n 2 

-2(1963 + 16^ 2 - 1552^ 3 )n + 64(29 - 2l£ 3 )] 
- 4C F C A [(4 + 45^ 2 - 48^ 3 )n 4 - 6(29 + 37^ 2 - 70^ 3 )n 3 + 3(228 + 139^ 2 - 408^ 3 )n 2 

-6(149 + 62^ 2 -242^ 3 )n + 4(95 + 33^ 2 - 150^ 3 )]} . (4.9) 

The corresponding (and more complicated) coefficients for the gluon case are 
C lT = ^^« / {cr 1 (^(77-15^ 2 -4(5 + ^ + ^(17-3C 2 )) 

- SC F 5 n2 - 2(4C F p - 3 \C 2 F + 1 6C& 2 ) 8„ 3 + (in 3 -^n 2 -2n + ^ Q n2 
-C n A 2 C F [16(5 - ^ 2 )n 2 - 8(29 - 5^ 2 )n + 8(20 - 3£ 2 )] 0„ 3 

jabc j 

- C A "- 3 a -^- fl 8 n (32n 2 - 96n + 64) ( 1 1 + 2^ 2 - 1 2^ 3 ) 0„ 3 

fta 

^n-3ol( 1 4 29 3 41 2 49 1 \ 

+ C; 3 p o 2 (-n 4 --n 3 + -n 2 --n + -)0„3 

+ C;- 1 [2(9 - 8^ 2 )n - 16(2 - £ 2 )] Q n4 - C^fo ( \n 2 + \n-l) %A 

+ C^Cj [20n 2 - 8(7 + ^ 2 )n + 2( 17 + 4£ 2 )] 9„ 4 

-c;- 3 c F (3 (n 3 -^ 2 + ^-i)e„ 4 - £c/c A ^- 2 po(^ 2 -^+£)0„4 

+ EC/C;-^ 1 [4(1 -2^ 2 )n + 6£-2(9-4^ 2 )]0„ 4 ) , (4.10) 

1=3 > 

rff = -j^-n f [l6C F ^n2 - 8(3C F (3 -12C| + 4CK2)5„3 

+ Q"- 1 ( ^(77 - 15^ 2 )n 3 - 8(16 - 3^ 2 )n 2 + |(61 - 21^)n - 8(7 - 3£ 2 )) 
+ C;- 2 (3 (4n 4 - i^n 3 - 24« 2 + - 36) 0„ 2 

- C n A 2 C F [16(5 - ^)n 3 - 8(39 - 7^ 2 )n 2 + 24(16 - 3£ 2 )« - 32(5 - £ 2 )] 6„ 3 
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n-2 

—1 

~3 

and finally 



r ]abc ] 

- C A "- 3 n abc fl 8 n (32n 3 - mn 2 + 160n - 64) (1 1 + 2£ 2 - 12^)6„3 

^n-3o2f 1 5 31 4 59 3 319 2 179 5\ 

+ C n F - 1 [2(9 - ^ 2 )n 2 - 2(37 - 24^ 2 )« + 4(19- 8£ 2 )] 6„ 4 
- C A "- 3 C F (3 (n 4 - ^ 3 + ^n 2 + 4n - 6) n4 

+ C A "- 3 C| [20n 3 -4(19 + 2^ 2 )n 2 + 2(43 + 12^ 2 )n - 2(17 + 8£ 2 )] 9 n4 

-c;- 2 po(^ 3 + ^ 2 -6n+8)e„ 4 - £c/c A "- / - 2 (3o(^ 3 -^ 2 +^-4£)e„4 

+ |4cf w [4(l-2^ 2 )n 2 + (6£-2(13-12^ 2 ))n-18£ + 2(21-8C 2 )]0„ 4 } (4.11) 



cif = ^,4cr 2 Po(2^-^-^ 3 + iii„ 2 -f „ + f )e„ 2 

-4(3C F po-10C 2 )5„ 3 + C;- 1 (^(77-15^)n 4 -^(127-27^ 2 )n 3 

+2(91 - 29^ 2 )n 2 - |(49 - 27^ 2 )« + y (5 - 3£ 2 )) 
-C A "- 2 C F [8(5 - ^ 2 )n 4 - 4(59 - 1 l£ 2 )« 3 +4(125 - 24^)« 2 
-4(114-25^ 2 )n + 4O(4-^ 2 )]0„3 

jabc j 

_ c n-3 abc fl 8^ l6n 4 _ ^ + 2QSn 2 _ 192n + 64 ) ( ! ! + 2 ^ 2 _ 12^ 3 )e„ 3 

^_ 3o2 / 1 6 13 5 79 4 1217 3 65 2 473 11 \ 

+ C7 3n — n n n H n 3 n 2 H n 0„ 3 

H0 V32 288 96 288 8 72 6 / 

+ C;- 1 [(9 - 8^ 2 )n 3 - (67 - 48^ 2 )n 2 + 4(41 - 22^ 2 )n -12(11- 4£ 2 )] 0„ 4 
+ C A "- 3 C| [lOn 4 - 2(29 + 2^ 2 )n 3 + 3(39 + 8^ 2 )n 2 -11(9 + A^ 2 )n + 8(5 + 3^ 2 ) 
- C;- 3 C F (3 {^-n 5 - 5 -n A - 3n 3 + ^-n 2 -19n+ 15) 6 n4 

-Cr 2 p (V-yn 2 + yn-18)e„4 

" E C/c;-^ 2 (3 ( V - in 3 + I (2£ - 5)n 2 - \ (3t - 1 )n + 9i) 6„ 4 



9;j4 



=2 
n-2 

+ £ C/c;-^ 1 [2(1 - 2^)n 3 + (3£ - 19 + 24^ 2 )n 2 - (2U - 63 + 44^ 2 )n 

^=3 

+36£-24(3-^ 2 )]0„ 4 l • (4.12) 
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5 NNLL resummation of the off-diagonal splitting functions 



Together with the mass-factorization relations of Section 2, the results of the previous two sec- 
tions facilitate the iterative determination of the respective coefficients D^- 2 ' and D^f- 2 ^ in 
Eqs. (11.31) . Eqs. (|1.5I) and Eqs. (11.61) to any order a". In this section we present the resulting expres- 

(n £<2) 

sions for the splitting functions in Mellin-/V space, from which the jc-space coefficients Dq g ' - 

and Dgq^- 2 ^ can readily be obtained by inverting the second relation in (12.21) . The corresponding 
results for the coefficient functions are discussed in the next section. 

The LL and NLL contributions to P qg and P gq can be expressed in a closed form in terms of 
a new class of functions with Taylor expansions in terms of Bernoulli numbers. Extending the 
definitions of Ref. lfT5l to k = 2, 3, . . . , these functions are given by 

B n are the Bernoulli numbers in the standard normalization of Ref. 11521 : i?2n+i = for n > 1 and 

Bn = 1 , B i = — , Bo = — , Ba = , Bf. = — , .... B 1 9 = , ... . (5.2) 

u ,1 2 , i 6 , 4 30' 6 42' ' 2730 

The functions (B±k for k > are related to <Bq, which can also be written as 
by 

—r {x k <B k ) = «0 , -rr $0 = -r ®-k ■ (5-4) 
ax K ax K x K 

Due to C,2n — > 1 for n — > °« the series (15.11) converge absolutely for all values on x. 

The numerical behaviour of the functions with k = 0, ±1,— 2 which enter our results 

below is illustrated in Fig. [U Similar to (Bq(x) - where this oscillation has been shown to continue, 
albeit in a much more irregular fashion, to much larger values of x 11531 - r Bk>o(x) oscillates around 
y = for positive x and y = —x/(k+l) \ for negative x. On the other hand, &k<o( x ) a l so oscillates 
around y = for positive x but around y = — x for negative x. As can be seen from the figures, the 
amplitude of these oscillations increases very rapidly with decreasing k. 

The resummed gluon-quark splitting function at large Af can now be written as 

ArP qg (Ar,a s ) = 2a s n f <B (a s ) 

+ a s 2 ln/YnJ(12C F -2p ) ^®-l(«s) + ^<B-i{a,) + (6C F - p ) ®l(«s) 
■ J I a s a s 

+ NNLL contributions + ... . (5.5) 
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Here and below we use (recall InN = lniV + y e ) 

a s = 4a s (C A - C F ) \n 2 N . (5.6) 
The corresponding result for the quark-gluon splitting reads 
NP gq (N,a s ) = 2a s C F $ (-a s ) 



+ ahnNCf 



;i2C F -6p ) ^®-i(-a s ) - ^-B_ 2 (-a s ) 



+ (14C F -8C A -p ) «i(-«s 
+ NNLL contributions + . . . . (5.7) 

The first lines of Eqs. (15.51) and (15.71 ) are the respective LL results, derived in Ref. [fT5l from 
relations equivalent to Eq. (13.121) . Unlike for the NLL corrections, we have not been able to find 
closed relations for all colour factors contributing to the NNLL terms. We therefore provide these 
results in the form of tables to order a s 18 which can be found in Appendix A. 

Our results for the fourth-order (N 3 LO) splitting functions agree with the predictions of Ref. 
ifTTl derived from the conjectured single-logarithmic large-jc enhancement of the physical evolution 
kernels for the system (F2, F^) of flavour-singlet structure functions. Furthermore they show the 
expected extension of the colour-factor pattern seen in those results to all orders in a s : the LL terms 
for both Pqg and P gq [recall Eq. (|l.ll) l are proportional to (C> — C A ) n , the NLL terms include 
at most two colour factors other that (Cf — Ca), and the NNLL terms of Appendix A involve 
(C F -C A )"- 2 or higher powers of (Cf — Ca)- It appears that generally all double-logarithmic 
contributions, \n k N with n + 1 < k < 2n, vanish for Cp = Ca which is part of the colour-factor 
choice leading to an 5\£ = 1 supersymmetric theory. 

It is clear already from the discussion of the fourth-order results in Ref. lfT71 (Figs. 12 and 
13) that the highest three logarithms of Pqg^ and Pgq^ can provide no more than a very rough 
indication of the size of the corrections beyond NNLO (note that the expansion in Ref. fT71 is in 
terms of InN instead of InN). Nevertheless it is useful to take a brief look at the numerical size 
and convergence of the resummation corrections to these quantities. This is done in Figs. 2 and 3 
at the reference point 

cc s (<2 2 ) = 0.2 , n f = 4 (5.8) 

used before, e.g., in Refs. [|8]4T2l. Depending on the precise value of a s at the Z-boson mass, this 
choice corresponds to a scale Q 2 ps 25 ... 50 GeV 2 typical for measurements of DIS HI. 

The figures show that the resummation corrections to the NNLO splitting functions are domi- 
nated by the third (NNL) logarithms which completely overwhelm the LL and NLL contributions 
(except, of course, for huge but practically irrelevant values of N). The relative corrections to the 
already small large-TV off-diagonal splitting functions are rather small, amounting to less than 2% 
and about 3% at TV = 20 for for P qg and P gq respectively. At least at the present level of logarithmic 
accuracy contributions beyond order (N 3 LO) are negligible. 



20 



« 1 1 1 1 1 1 

\ 


i i i i i 

np^n) : 

- 


\ 




\v 








■ NNLO 




1 LL 




• NLL 




■ NNLL 

i 


1 



10 



20 30 

N 



0.08 



0.075 



0.07 



0.065 



0.06 




Figure 2: The Mellin-/V splitting function P qg and P gq , multiplied by N for display purposes. Shown 
are the LL, NLL and NNLL large-Af resummation corrections to the complete NNLO results BUl. 
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Figure 3: The contributions of the various orders in oc s to the resummation corrections atN = 20. 
The LL, NLL and NNLL terms of order a" are added at the corresponding values of the abscissa. 
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6 NNLL resummation of the coefficient functions 



The resummed results for the coefficient functions C2, g and C^ q are more complicated than those 
for the splitting functions already at the leading-log level. Hence it is not surprising that we are 
not able to give closed NNLL expressions for these quantities either at this point. Postponing the 
NNLL contributions to Appendix B, the results are 
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Here the first lines are the respective LL results of Ref. |fT5ll . As in these terms, the (Ca — Cp) 
denominators are generally cancelled by corresponding numerator factors as can be seen by ex- 
panding all functions in powers of 0C S . Unlike for the splitting functions the double-logarithmic 
contributions to the coefficient functions do not vanish for Cp = Ca- However, they can be ex- 
pressed in terms of exponentials in this case since all <B -functions have the argument (|5.6I) . 
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Our results for the quark coefficient function for Fi completely agree with those of Ref. [|29ll . 
The complexity of the N"LL resummed expression for the gluon coefficient function C^ g is similar 
to that of C2, g at N" _1 LL level. Hence it can be written down in a closed form at NNLL accuracy, 



N z C L , g (N,a s ) = Sa s n f exp(2a s C A ln z N)+4a s C F NC2^(N,a s 
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+ ~ a 2 ln 4 /V p c A ((3 + 16C A - 4C F ) + - a] ln 6 ^V(3 2 C A 2 

Here the first term is the LL result also obtained already in Ref. [|29l . The next term, where 
Cfjr stands for the first line of Eq. (16.11) , and the second line form the NLL contribution. The 
remainder of Eq. (16.31) collects the NNL logarithms, where C^g L represents the second to sixth 
line of Eq. (16.11) . Furthermore this is the only NNLL contribution including the charge factor 

flf x = (e) 2 /(e 2 ) with (e k ) = nf 1 l,Zi e i ( 6 - 4 ) 

(where et is the charge of the z'-th effectively massless flavour in units of the proton charge) arising 
from diagrams with the colour factor d abc d a }, c /n a = 5/48 nj in QCD where the two (neutral) gauge 
bosons couple to different quark loops of the gauge-bosons gluon forward amplitude f^l fTOlfTTTl . 

The NLL and NNLL contributions to Eq. (16.31) for C F = and flf x = agree with the previous 
all-order result (H in this gluonic 'non-singlet' limit to which C^ LL and hence the special func- 
tions (15.11) do not contribute. The complete fourth-order coefficient function (for W-exchange, i.e., 
without the flf^ = part) has been predicted in Ref. [54] from the physical evolution kernel for the 
system (F2, Fi) of flavour-singlet structure functions together with the four-loop splitting-function 
results of Ref. [fTTl . The present results are in full agreement also with that prediction. 

Eqs. (15.51) and (15.71) for the resummed splitting functions together with their counterparts (16.11) - 
(16.31) for the coefficient functions and the corresponding simpler results for C2. q , C* g and Cl a 
facilitate an assumption-free NNLL calculation of the physical kernels for (i^i 7 ^) and (F2,F F ) to 
any order in oc s . It turns out that their highest double logarithms, as far as they can be determined 
from available fixed-order results now, do indeed vanish. Hence the conjectures of Refs. lfP71l28l 
54l are proven by our present calculations for the leading TV -1 large-TV contributions. 



The numerical size and a s -convergence of the LL, NLL and NNLL resummation corrections 
to the respective third-order results are illustrated in Fig. 4 for C2 ;g and Figs. 5 and 6 for the 
coefficient functions of F F . For brevity we do not show the corresponding results for C^q which 
is of theoretical but not phenomenological interest. The corrections are dominated by the NNLL 
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Figure 4: Left: the Af-dependence of the third-order (N 3 LO) and resummed gluon coefficient 
function for F 2 at the reference point (15 . 81 ) . Right: the contributions of the fourth to the tenth 
orders in oc s , added at the corresponding values of the abscissa, to those results at Af = 20. 



contributions, suggesting that they underestimate the impact of at least very high orders. Already 
the known terms, e.g., at N = 20 for the reference point (15.81) . are sizeable for Ct,q with about 15%, 
large for C2, g (about 35%) and huge for C^g with about 100%. The gluonic quantities receive 
significant contributions from the fourth to sixth order in oc s . With values of 10.2 and 49.5 the 
forth-order coefficients are comparable to the corresponding (averaged) Pade estimates, 

- c 2 , g (N = 20) = 0.127a s + 0.642 a s 2 + 2.76 a s 3 + 12p ad6 a s 4 + . . . , (6.5) 
NC L , g (N = 20) = O.llOcts + 1.240a s 2 + 9.51a 3 + 65 Pad6 a s 4 + ... . (6.6) 

It is obvious from Eqs. (13.181 ) - (13.211) that our method can be applied as well to non-leading 
terms in the expansion in powers of ( 1 — x) . This is especially interesting for the subleading ( 1 — x) 
or 1/N contributions to diagonal coefficient function such as C2, q and C^g. Due to the stable 
form (11.21) of the diagonal splitting functions, the corresponding N 3 LO corrections to the non- 
singlet structure functions are known except for the single- and non-logarithmic contributions. 
Consequently, as discussed below Eq. (13.211) , we are able to resum highest four A^ _1 logarithms. 

The resummation of the Af -1 non-singlet coefficient functions for ^1,2,3 has been inferred in 
Ref. [28J from the behaviour of the physical evolution kernels, obtaining complete NNLL results 
and the N 3 LL corrections up to one undetermined number called ^dis 4 - Hence it was not necessary 
to perform a very cumbersome all-order N 3 LL calculation analogous to Section 3. Instead we have 
verified those previous results to order aj and determined the hitherto missing parameter, 

^ DI§4 = 100/3 in Eqs. (5.24) - (5.26) of Ref. . (6.7) 
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Figure 5: The Mellin-Af space quark (left) and gluon (right, multiplied by AO coefficient func- 
tions for Fl at the reference point (15.81) . Shown are the cumulative LL, NLL and NNLL large-Af 
resummation corrections to the third-order (NNLO) results. 
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Figure 6: As Fig. 3, but for the coefficient functions Ci, q and C^ g shown in Fig. 5 above. 
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7 Summary and outlook 



Unlike the case of the dominant (1— x)T /N large-*/ large-Af terms, only a small amount of 
(published) research has been devoted until recently to the all-order structure of 1 /iV-suppressed 
threshold contributions to the coefficient functions for the structure functions in deep-inelastic scat- 
tering (DIS) and crossing-related (semi-) inclusive quantities in perturbative QCD. See Ref. 11551 
for a well-known leading-log conjecture, and Refs. Il56l [571l for early studies of the longitudinal 
structure function F^. The study of these contributions is however not only theoretically interest- 
ing but also phenomenologically relevant, e.g., for assessing the kinematic region, different for 
different processes, see Ref. fl2l . in which the N° terms and their soft-gluon exponentiation can 
be used as a quantitative substitute for the full coefficient functions. 

Consequently several groups have addressed this issue with various approaches in the past few 
years, for work by others see Refs. |j58]46TTl and ll62l - [64ll . However, to the best of our knowledge, 
this research has not yet led to explicit all-order predictions for complete and exact coefficients 
of the next-to-leading and next-to-next-to-leading logarithms, or even the leading logarithms to 
quantities such as the gluon contribution to most important structure function F%. 

In Refs. Il28~ll2~9l a resummation of the highest three l/Af-suppressed logarithms (actually the 
l/N behaviour is more formal than real for practically relevant values of N) for various quark 
coefficient functions has been obtained by studying the large-Af behaviour of non-singlet physical 
evolution kernels. These quantities express the scaling violations of a physical quantity in terms of 
the same physical quantity and hence do not depend on the scheme for the factorization of the mass 
singularities. It turns out that the non-singlet physical kernels show an only single-logarithmic 
large-Af enhancement also beyond the dominant N° contributions (where this feature is a simple 
consequence of the soft-gluon exponentiation, see Ref. Il65l ) up to at least the next-to-next-to- 
leading or next-to-next-to-next-to-leading orders (NNLO or N 3 LO). The results of Refs. 112811291 
are based on the conjecture that this behaviour continues to all orders in the strong coupling oc s . 

Completely analogous observations were made lfT7ll54ll for the physical-kernel matrices of the 
systems (F2,F§) and (F2,F£) of flavour- singlet structure functions [66-68]. Unlike in the non- 
singlet case, recall Eq. (11.21) . the singlet kernels receive double-logarithmic contributions from 
both the splitting functions (11.31) and the coefficient functions in Eqs. (11.51) and (11.61) . Conse- 
quently the single-logarithmic enhancement of the physical kernels can only provide one all-order 
constraint between two quantities. Definite fourth-order predictions of the highest three logarithms 
are possible though: first for the N 3 LO splitting functions [fTVTl using the diagram calculations of 
the corresponding three-loop coefficient functions [fTQllTTllTvTl . and then, using these predictions, 
for the N 3 LO (fourth-order) coefficient functions for Fl ll54l . 

While we expect the results of Refs. H17[|28[|29ll54l to be the final word on the predicted coeffi- 
cients, clearly more work is required to put them on a firmer theoretical footing and to extend them 
to all orders also for the off-diagonal splitting functions and flavour- singlet coefficient functions. 
At the leading-logarithmic level this was done in Ref. |[T5l by finding the now expected all-order 
iterative structure of the unfactorized partonic structure functions (forward-Compton amplitudes). 
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In the present paper we have extended these results to the next-to-leading and next-to-next-to- 
leading (NLL and NNLL) logarithms. Note that this is the counting of a resummation, not that 
of a more powerful exponentiation [fT8]422| : the present predictive power in terms of higher-order 
coefficients corresponds to that of a next-to-leading logarithmic exponentiation, see, e.g., Ref. Il2~5ll . 

In order to achieve this, we have employed two independent methods. The first is a direct 
generalization of the amplitude iteration of Ref. lfT5l to higher logarithmic accuracy, presented at 
NLL accuracy in Eqs. (13.141) - (13.161) above. The second, worked out to NNLL accuracy in the rest 
of Section 3 and Section 4, is conceptually extremely simple and revealing but not iterative, since 
all coefficients are determined order-by-order in a s without any explicit reference to complete 
lower-order amplitudes. This method is based on only the form of the D-dimensional phase space 
for inclusive DIS and the all-order mass -factorization formula (guaranteed by the operator-product 
expansion) for the unfactorized structure functions. In this way we have been able to verify all DIS 
predictions in Refs. H17[|28[|2"9ll54"l . extend them to the fourth logarithms for F2, ns and to all orders 
in oc s for the highest three logarithms of the off-diagonal splitting functions and the coefficient 
functions C2, g , C^ g and C^ q. Our results prove the non-singlet conjecture of Refs. If28ll29l and 
show that also the above singlet physical evolution kernels are single-log enhanced to all orders. 

The resulting all-order off-diagonal splitting functions and the above coefficient functions can 
be written down in a closed all-order form to NLL (NNLL for C^.g) accuracy in terms of the 
apparently new special functions introduced in Ref. [15J. We did not find (so far) similar relations 
for all colour factors of P qg and P gq and for C2. g and C^.q at the NNLL level, but can analytically 
determine their coefficients to 'any' order in a s . The splitting-function results show a particular 
colour structure with all double logarithmic contributions, a" In (1— x) with I > n, vanishing in 
the supersymmetric limit C> = Ca- This pattern does not hold for the more complicated coefficient 
functions which however can be expressed in terms of exponentials in the limit. 

Numerically the resummation corrections to the a 3 splitting functions appear to be rather small 
and quickly converging: at the present level of accuracy only the fourth-order terms have an im- 
pact. The corresponding corrections for C^ q are larger and receive relevant fourth- and fifth-order 
contributions, but they are still small compared to those up to order a s 6 for C2 ;g and especially C^g. 
Taking into account that the highest three logarithms most likely underestimate the corrections at 
high orders of a s , these results reinforce the third-order findings of Ref. [1 1 J which indicated that 
the perturbative expansion of Ci, g is not well-behaved even at moderately large N in the important 
region of scales Q 2 ~ 25 ... 50 GeV 2 corresponding to 0C S (<2 2 ) ~ 0.2. 

Beyond first estimates of the numerical impact of higher orders at large N, our results will 
prove very useful once the fixed-moment calculations of Refs. [|5]-[7] have been extended to the 
fourth order in a s . First results have already been presented on four-loop second moments and 
sum rules Il2"6ll69ll . hence this extension can be expected in the foreseeable future. The main way 
to use a couple of moments of, e.g., the N 3 LO splitting functions will be by effective x-space 
parametrizations analogous to those of Refs. H70U71H at NNLO. Already knowing the coefficient of 
\n e (l—x) for £ = 4, 5, 6 will then considerably assist achieving a decent accuracy from a limited 
number of moments. The situation is analogous for the coefficient functions. 
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There appears to be no reason why our determination of the highest double logarithms cannot 
be extended beyond the N~ l (N~ 2 for Ct.g) terms. An all-order calculation would obviously be 
very cumbersome given the form of some of the intermediate expression in the present limit, but 
a fourth-order calculation would definitely be feasible. We would expect to recover the results 
of Refs. 1(1711281 for the coefficient of the highest three logarithms at all-powers of (1 — x) of c^q 
with a = 2, 3, L, Pqg and P g q in this manner, and to derive corresponding new results for ci 4 g 

(4) 

(a = 2, L) and cl '. The main application of such results would be as a check of an all-x fourth- 
order calculation of the splitting functions and DIS coefficient functions which, however, we do 
not expect for the near future. 

Due to the similar phase-space integrations J72l|73l, the present approach can be carried over 
directly to the case of semi-inclusive e + e~ annihilation [1J (but not, unfortunately, to the Drell- 
Yan process and inclusive Higgs production in pp collisions). Only the diagonal NNLO 'timelike' 
splitting functions for parton fragmentation have been determined up to now [174117511 . hence a 
NNLL resummation is not yet possible in this case as the third logarithms of the NNLO splitting 
functions are a necessary input. Therefore this case will be addressed in a future publication after 
the determination of the off-diagonal timelike splitting functions 11761 . 

Let us finally stress again that the present double-logarithmic resummation is not relying on 
any specific large-Af structure beyond the general form of the phase space integrations. One may 
hope that, as in the case of the soft-gluon exponentiation of the N° coefficient functions, such 
additional structures can be found, e.g., using the new approach of Refs. [|59l[60ll or by improved 
applications of the soft-collinear effective theory to DIS, see, e.g., Refs. Il77148~0l . It should then 
become possible to resum also lower logarithms analogous to the standard threshold resummation. 
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Appendix A: Large-iV splitting functions at NNLL accuracy 



Here we provide the expansion coefficients of the highest three logarithms of the off-diagonal 
splitting functions, defined via 
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with Pqg(N) = Pqg/np Pgq(N) = Pgq/CF and Caf = Ca~Cp. It is understood that the coefficients 
on the r.h.s. depend on the splitting function under consideration. Dnll.i vanishes for P qg , but not 
for P gq . The coefficients of the leading and next-to-leading logarithms, for which closed expres- 
sions have been written down in Eqs. (15.51) and (|5.7I) . are given up to order a s 18 in Table 1. Their 
counterparts for the third logarithms are provided in Table 2 for P qg and Table 3 for P gq . The closed 
expression leading to most of these coefficients is not known yet. 



Appendix B: Large-N coefficient functions at NNLL accuracy 

Finally we turn to the NNLL contributions to the off-diagonal coefficient functions C2, g and C^q. 
The more complicated colour structure, with 3n — 3 contributions at order n, precludes a represen- 
tation in the form of Tables 2 and 3. Instead we provide the resulting expression at order ocf and 
0C S 5 for general SU(N) colour factors, and then the six-figure values for QCD to order cc s 12 which 
are more than sufficient for numerical applications. In all cases we include also the LL and NLL 
terms, for which closed expression have been given in Eqs. (16.11) and Eqs. (16.21) above. 

The fourth- and fifth-order contributions to the gluon coefficient function for F% are given by 
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Table 1 : The coefficients of the leading and next-to-leading large-Af logarithms of the off-diagonal 
splitting functions P gq and P qg as defined in Eq. (IA.1I) up to the 18-th order a s = oc s / (4tt). Note 
the respective appearance of the numerators 691 and 3617 from n = 12 and n = 16 which clearly 
signals the presence of the Bernoulli numbers (15.21) also in the NLL coefficients. 
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36728 
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91 
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12 
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733288 


71989 
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3617 
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Table 2: The coefficients of the next-to-next-to-leading large-Af logarithms of the off-diagonal 
splitting function P qg , as defined in Eq. (|A.1I) . up to the 18-th order a s = oc s / (4tc). 
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80 
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67 


9 
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4~ 


T 


35 
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~y 
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Table 3: As Table 2, but for the splitting function P m . 
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The corresponding results for the quark coefficient function for F§ read 
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Our notation for the Tables 4-6, where we include for the convenience of the reader, is 

^5 (ao = 
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with C2 ig = C 2 , g /n, and C^.q = C^^/Cp. fl 8 n has been defined in Eq. (16.41) . 



(B.6) 



Finally it is worthwhile to note that a closed (if presumably rather lengthy) expression for the 
NNLL contributions to C2 ig and CVq can be derived once such an expression has been obtained 
for the corresponding contributions to the splitting functions in Appendix A. At that point all 
quantities but C2, g and C^q entering the vanishing off-diagonal NNLL elements of the physical 
evolution kernel 



dF f . dC 
dh^ = ^ {as) Ta, 



($(a s )—C- l +CPC- l )F = KF 
V da* J 







NNLL 



with the standard matrix P of the singlet splitting functions and 



F 2 
Fa, 



C 



C 2,q C 2,g 

H>,q H>,g 



K 



K22 
K§2 



(B.7) 



(B.8) 



will be known, and Eqs. (|B.7D can be solved for the NNLL parts of C2, g and q . In fact, we have 
applied the analogous NLL procedure to find the closed expression given in Eqs. (16.11) and (16.21) . 
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3 


11.9259 


103.944 


1.61728 


285.481 


13.7942 





4 


16.7874 


248.091 


5.73937 


1493.50 


82.2181 


0.489712 


5 


19.5455 


419.726 


11.6100 


4059.62 


248.604 


2.43393 


6 


19.3641 


561.667 


17.4564 


7676.98 


508.882 


6.25137 


7 


16.5861 


624.621 


21.2105 


11287.1 


797.203 


11.2984 


8 


12.4589 


592.811 


21.5942 


13619.9 


1014.87 


15.9871 


9 


8.31829 


489.124 


18.8489 


13915.5 


1085 50 


18.5769 


10 


4.99558 


356.178 


14.3660 


12296.6 


997.538 


18.2401 


11 


2.72610 


231.848 


9.70788 


9552.90 


801.374 


15.4587 


12 


1.36329 


136.369 


5.89075 


6612.08 


570.865 


11.5051 



Table 4: The LL, NLL and NNLL coefficients of C 2 , g in QCD, as denned in Eq. (|B3T) . to the 12-th 
order in a s = oc s / (47l) . All these coefficients further decrease at higher orders and tend to zero in 
the infinite-order limit. 
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5 


123.863 


1261.66 
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11928.5 
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25.3455 


6 


151.299 


2363.55 


148.156 


27287.3 


2950.46 


68.9743 


7 


154.905 


3455.21 


215.475 


49475.2 
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139.478 


8 


136.254 


4124.46 


256.338 


73284.1 


8381.24 


222.015 


9 


104.882 


4147.64 


257.174 
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10551.2 


289.008 


10 


71.6724 


3594.93 


222.527 


95974.4 


11279.5 


316.347 


11 


43.9873 


2732.85 


168.950 
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10424.5 


297.446 


12 


24.4770 
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114.098 


70891.3 


8453.50 


244.331 



Table 5: As Table 4, but for the coefficient function C^.q. 
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28.4444 
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10 
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620424. 
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266.24 


691.2 


11 


133.303 
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529638. 


38847.9 


199.68 


658.286 


12 


72.7106 


7717.03 


296.229 


398936. 


30392.6 


133.12 
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Table 6: The LL, NLL and NNLL coefficients of C L , g in QCD, as denned in Eq. (lEToT) . to the 12-th 
order in a s = oc s / (4rc). Also these coefficients further decrease at higher orders and tend to zero in 
the infinite-order limit. The fifth column represents the fl g n term absent in charged-current DIS. 
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